Reconstruction of dynamical equations for traffic flow 



S. Kriso^, R. Fricdrich^ , J. Peinke^ and P. Wagner'^ 
^ Fachbereich 8 Physik, Universitat Oldenburg, 26111 Oldenburg, Germany 
^ Institut fiir Theoretische Physik und Synergetik, Universitat Stuttgart, 70550 Stuttgart, Germany 
^ Deutsches Zentrum fiir Luft- und Raumfahrt (DLR) e.V., Institut fiir Verkehrsforschung, Rutherfordstrasse 2, 12489 Berlin, 

Germany 
(February 2, 2008) 



Traffic flow data coflected by an induction ioop detector on 
the fiigfiway close to Koin-Nord are investigated with respect 
to their dynamics including the stochastic content, fn partic- 
ular we present a new method, with which the flow dynamics 
can be extracted directly from the measured data. As a result 
a Langevin equation for the traffic flow is obtained. From the 
deterministic part of the flow dynamics, stable flxed points 
are extracted and set into relation with common features of 
the fundamental diagram. 

I. INTRODUCTION 

By increasing the number of licensed vehicles on our 
roads it becomes more and more necessary to reduce ap- 
pearances of traffic congestion. Getting higher capacities 
of highways means looking for optimized flow rates of 
cars. In order to do so, it is necessary to investigate the 
complex system of traffic flow and to understand its reg- 
ularities. In a subsequent step, intelligent traffic control 
systems may use these laws to influence the traffic flow 
and thus to increase the highway's capacity utilization. 

For our investigation a large amount of traffic flow data 
were collected at the highway near Koln-Nord (Germany) 
over more than one week. For each car crossing an in- 
duction loop detector the following data were recorded: 
(i) the time, when the car crossed the detector, with an 
resolution of 1 sec, (ii) the type of car (passenger car, 
truck), (iii) the lane number, (iv) the velocity, truncated 
to an 8-bit integer (0...255 km/h), accuracy of approx- 
imately 3 %, (v) the length of the car in meters (8-bit 
integer, not calibrated) and (vi) the distance to the car 
driving ahead (an integer in the range 0. . . 999 m, not 
calibrated). The lanes are labeled from A to C, where 
lane A is the right driving lane where mostly trucks are 
found (because of a German law) , and lanes B and C are 
the fast drivers' lanes. 

In the following we present a new method to derive 
from the measured data dynamical equations for the traf- 
fic situation. We start with the common presentation of 
the fundamental diagram. Next, we evaluate the deter- 
ministic and the stochastic content of of the traffic dy- 
namics by means of a Langevin equation. Finally, we give 
an interpretation of the determined Langevin equation. 



II. FUNDAMENTAL DIAGRAMS 

Here and in the following we calculate the car density 
k [km^^] and the flux state q [h~^] of each car using the 
measured data: the velocity v [km/h], the length of the 
car I [m] and the distance d [m] to the car ahead: 



q = k-v . (2) 

To avoid an overload of the presentation we restrict 
the diagrams and the following calculations to two cases. 
Firstly, only the traffic of a single lane C, secondly the 
cumulative traffic of all three lanes A, B and C are con- 
sidered. Furthermore we calculate for each density state 
k the mean flux state 

(^)^ = ;^E'^w- (3) 

The obtained fundamental diagrams (qi;ki) and 
(<q>;k) are shown in Fig. |^. We note that for this pre- 
sentation no significant difference in the traffic dynamics 
of one lane and the cumulative dynamics of three lanes 
can be detected. In both diagrams we find a maximum 
flux in free traffic flow of g,nax = 8000 cars/h, according 
to the results of For the flux out of traffic jams we 
find in both cases gout = 5000 cars/h. So we have a ratio 
of 77 := 'Zmax/'i'out = 1.6, which meets the value 77 « 1.5 
found in quite well. 

III. LANGEVIN EQUATION FOR THE TRAFFIC 
FLOW 

In order to grasp the underlying dynamics of the traf- 
fic flow (one lane and all three lanes), we utilize a new 
method to analyze the traffic data more extensively. In 
particular, the iterative dynamics of the traffic state Xjy, 
given by the velocity v^y and the flux of the iV-th car, 
as a function of traffic state X7v_i of the N — 1-th car is 
investigated. Note, that other state variables could have 
been choosen as well. 

For the iterative dynamics of the traffic state xjv 
(which may be taken in a generalized case as a r- 
dimensional variable) we propose a description by a sta- 
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tionary Langevin equation, taking into account a com- 
bination of deterministic and random (noisy) forces cf. 



where the indices i and j denote components of the mul- 
tidimensional variables, N the time, and F^jv are r x N 
independent Gaussian noise variables with zero mean and 
with variance 2, i.e., 



(5) 



The central part of our following work is that it is pos- 
sible to determine the functions hi and gij directly from 
empirical data. Taking as the Ito presentation of the 
stochastic process the following relation to the Kramers- 
Moyal coefhcients can be given. 



(6) 
(7) 



fc=i 



where D^^^ and D^^^ are called drift and diffusion coeffi- 
cient. These coefficients can be evaluated by the condi- 
tional moments 



■'(x) = {Xi^N+l - Xi) 

Dlj'{x) = {{xi^N+i ~ Xi){xj^N+i - Xj)) 



(8) 
(9) 



Recently it has been shown that with the analogous 
definition of these Kramers-Moyal coefhcients it is pos- 
sible to reconstruct from time continuous dynamics the 
underling stochastic differential equation [||. 

Before presenting our results on the dynamics we want 
to comment on the validity of this ansatz to describe the 
traffic flow by the Langevin equation (|^). This ansatz 
implies that the dynamics is in the class of Markovian 
processes, i.e. the system does not have a memory. This 
can be tested by conditional probabilities 



p(xn|xn-i, . . • ,XN-m) = p(xn|xn-i) 



(10) 



or by the necessary condition of the Chapman- 
Kolmogorov equation 

p(xN|xN-r) = ^ p(xN|xN-s)p(xN-s|xN-r) , (H) 



where r > s. From our data, conditional probabilities 
have been evaluated and the validity of the Chapman- 
Kolmogorov equation was found for the iterative dynam- 
ics of both quantities, the velocity and the flux. If this 
Markovian property holds, the inherent noise of the dy- 
namics (F) can be taken as J-correlated. It should be 



noted, that even in the case where the noise is not 6- 
correlated, the deterministic part of the dynamics can 
be reconstructed from given data (^, as we found by 
analysing numerically generated test data, . 

As expressed by (^, the knowledge of the conditional 
probabilities p(xn|xn_i) provides the basis to estimate 
the Kramers-Moyal coefficients from the traffic data. 
First we consider the simplified case of the onedimen- 
sional dynamics of the velocity only. The results for and 
D'^'^^v) for the traffic of one lane and for the cumulative 
trafhc of all three lanes are shown in Fig. ^j. The one di- 
mensional deterministic dynamics can also be expressed 
by the potential $d(v), defined as - ^ = D(i)(v). The 
corresponding potentials are shown in Fig. |^. From these 
results three noticable velocities Vi « 37 km/h, V2 « 75 
km/h and V3 « 107 km/h appear, which allows to iden- 
tify the so called congested flow for v < vi, correlated 
flow for vi < V < V2 and the free flow for v > V2 [^,0, 
respectively. Note these velocities can be defined as fixed 
points (D^^) (v) — 0) of the deterministic part of the cu- 
mulative traffic dynamics (see Fig. 2b). 

For the traffic dynamics of lane C we find in the con- 
gested and in the correlated regime metastable trafhc 
states, the deterministic drift term Z?^^^ gets zero over 
finite intervals. This corresponds to the plateau struc- 
ture in the potential, see Fig. 3a. A clearly different 
behaviour is found for the cumulative traffic dynamics of 
all three lanes, see Figs. ^ and |^b. The different flow 
regimes are seperated by two fixed points at vi and V2- 
The slope of these fixed points defines the stability, thus 
the congested and the correlated flow regimes are sepa- 
rated by a stable fixed point, whereas the correlated and 
the free flow regime are separated by an instable flxed 
point. For the free flow, in both cases of one lane or 
three lane traffic a stable fixed point is found at V3, cor- 
respondingly the drift potential has its local minimum. 
Because of a speed limit of 100 km/h given on the in- 
spected highway we see a great increase of the potential 
for V > V3: the faster a car is driving, the stronger the 
attraction is to the potential's minimum |^ . 

To get an understanding of the real traffic dynamics 
grasped by these drift coefficients or drift potentials, the 
additional noise has to be taken into account. In Fig. ||c 
and d the corresponding magnitude of the noise are ex- 
pressed by the evalutaed diffusion coefficients D'-^\ The 
noise will now cause transitions between different flow 
states. For the traffics dynamics of one lane, the noise 
will effect larger fluctuations as it is the case for the cu- 
mulative traffic dynamics, which has two clear minima 
in the potential. A further interesting detail is that the 
magnitude of Z?^^^ has a minimum around the stable fixed 
point at V3. This indicates a pronounced stability of this 
trafhc situation. 

Next we present the results of a higher dimensional 
analysis by taking x with the components Xi = v and 
X2 = q- Now also the drift coefficient D^^^ becomes a 
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vector depending on v and g, as shown in Fig. ]^ These 
resuhs were obtained by binning the velocity flux data 
into a 25 X 25 matrix, corresponding to a binning of the 
velocity into intervals of Au = 5 km/h. 

The solid lines in Fig. ^ show the states where we have 
no drift of the flux component: Dq^' = 0. On these lines 
we find only a velocity drift with a constant flux. In ac- 
cordance to Fig. ||a we find in Fig. ||a for slow velocities 
mainly no velocity drift, in Fig. ^ (like in Fig. ||b) there 
seem to exist stable velocity drift states at the same ve- 
locities {vi ~ 37km/h,107 km/h). The topology of the 
instable fixed point gets now a saddle point which is at- 
tractive for larger and smaller flux values but instable in 
the direction of larger and smaller velocity values. Again 
clear differences of the dynamics of on lane and three 
lanes is found. 

IV. DISCUSSION AND CONCLUSION 

By the investigation of traffic flow data as an itera- 
tive stochastic process we were able to calculate from the 
given data the 1-dimensional drift and diffusion coeffi- 
cients and thus to find the deterministic and stochastic 
part of the corresponding Langevin equation. We are 
able to find stable, metastable and unstable states (fixed 
points) in the deterministic part of free, correlated and 
congested traffic flow. For a fully description of the whole 
dynamics also the diffusion coefficient has to be taken 
into account, which provides transition probabilities be- 
tween the different (meta-)stable states. Without this 
noisy part of the traffic flow dynamics, the stable states 
would never be left, i.e. a congestion would stay forever 
if once prepared. 

To see the dependency of velocity and flow from each 
other, the investigations were expanded to a higher di- 
mensional analysis. In this case we find the deterministic 
and stochastic part of the 2-dimensional Langevin equa- 
tion. Now we are able to identify stable velocity and fiux 
states of the deterministic part. 

Interestingly, the results found in this study are in 
agreement with empirical investigations that have iden- 
tified three phases of traffic |^ , together with transitions 
connecting theses phases. Especially the transition from 
free flow to correlated flow (or synchronized flow in the 
terminology of ^) has similarities. 

Finally we want to point out, that we presented here 
a new method to analyse traffic data with respect to a 
derivation of dynamical equations from pure data analy- 
sis. Furthermore we could show that our method provides 
more insight into the traffic dynamics than the presenta- 
tions of diagrams like the fundamental diagram. A clear 
difference in the dynamics of one lane and the dynam- 
ics of cumulative three lanes was found. Our analysis 
provides evidence of the presence of fixed points, which 



are of practical importance if a control of a traffic should 
be achieved. At last one may conclude that this method 
will be helpful to perform a more thorough comparison 
between traffic flow models and empirical data. 
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three lanes A,B, and C. The dots represent the measured traf- 
fic data, the solid line corresponds to the mean flux (g)fe. 
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FIG. 1. Fundamental diagram of (a) lane C and of (b) all 
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FIG. 3. Corresponding potentials for the deterministic dy- 
namics given by the drift coefficients in Fig. ^. 
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FIG. 2. Drift and diffusion coefficients for the traffic of one 
lane C, (a) and (c) and of all three lanes (b) and (d). 
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FIG. 4. Deterministic dynamics of the twodimensional 
traffic states (g, v) given by the drift vector for (a) one lane C 
and (b) all three lanes. Bold dots indicate stable fixed points 
and open dots saddle points, respectively. 



4 



